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Abstract 

This contribution covers the topics presented by the authors at the 
"Fundamental Problems of Turbulence, 50 Years after the Marseille Con- 
ference 1961" meeting that took place in Marseille in 2011. It focuses 
on some of the mathematical approaches to fluid dynamics and turbu- 
lence. This contribution does not pretend to cover or answer, as the 
reader may discover, the fundamental questions in turbulence, however, 
it aims toward presenting some of the most recent advances in attacking 
these questions using rigorous mathematical tools. Moreover, we consider 
that the proofs of the mathematical statements (concerning, for instance, 
finite time regularity, weak solutions and vanishing viscosity) may contain 
information as relevant, to the understanding of the underlying problem, 
as the statements themselves. 

MSC Subject Classifications: 35Q30, 35Q31, 76D03, 76D09. 
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limit. 

1 Introduction 

One of the main goals of this contribution is to underline the relationships be- 
tween, on the one hand, standard mathematical questions and, on the other 
hand, classical issues in turbulence theory. However, it is important to em- 
phasize that there are many different mathematical approaches aiming toward 
understanding the complex phenomenon of turbulence, and that only some of 
these aspects are stressed in this contribution (see also the other contributions 
in this special issue; and for a general user friendly introduction (for mathe- 
maticians and engineers) the book of Frisch j42j ) . 
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A nice survey of the dynamical systems approach to turbulence and statis- 
tical solutions of Navier-Stokes equations is presented in the monographs [27] . 
[4"T] . and in Foias [33] • Moreover, rigorous conditions for the main tenets of the 
Batchclor-Kraichnan-Lcith theory of 2D turbulence, namely the energy power 
law, enstrophy cascade, and dissipation law can be found in [3J, [35], [JU] and 
references therein. Furthermore, recent interesting surveys concerning analyti- 
cal studies of Euler equations can be found in [8J , [22] , [23] , [46] and references 
therein (see also [T5], [12], [1H] and references therein for recent surveys concern- 
ing computational studies of Euler equations). 

Starting from the macroscopic Navier-Stokes-Fourier system for compress- 
ible fluid flows, we will recall the derivation of the equations of motion for 
incompressible fluid flows, by taking the small Mach number limit. Then by 
another new scaling argument, which relies on the introduction of the Reynolds 
number, we will introduce the Euler equations as the formal zero viscosity limit 
of the Navier-Stokes equations. Historically, the Euler equations for ideal flu- 
ids were the first to be derived (Euler 1757), viscous effects being added only 
later to give the Navier-Stokes equations. These viscous effects are, as it will be 
stressed in the paper, essential for understanding the interaction between solid 
bodies and fluid flows. Moreover, recent mathematical results indicate that a 
more profound understanding of Euler equations seems to be essential for mak- 
ing any significant progress in this field. Therefore, section [3J is devoted to the 
mathematical theory of Euler equations with, on the one hand, emphasis on 
the classical local well-posedness of the problem in the class of smooth solutions 
and, on the other hand, on the more recent notions of weak solutions. We will 
then return, in section 4, to the Navier-Stokes equations, mostly to consider 
the boundary value problems and the behavior of the solutions, in particular 
boundary effects, when the viscosity tends to zero. 

2 Preliminaries 

Our approach follows three threads. First, the notion of the Cauchy problem for 
an underlying evolution equation, which means the study of the corresponding 
solutions of this equation, for t > to, given their state at time t = to (often to is 
taken equal to 0). Here the relevant issues are: the stability or instability of the 
solutions, their qualitative long-time behavior, and eventually the reliability of 
this equation for modeling the underlying phenomena at the relevant scales and 
under the appropriate conditions. 

The second thread is the notion of weak solutions. This idea preexisted in 
mechanics and physics in the context of, for instance, calculus of variations; 
however, in the previous century it became systematically used with the devel- 
opment of functional analysis and the introduction, by Gelfand, Schwartz and 
Sobolev, of the notion of distributions (which are also called, by some people, 
generalized functions). In a nutshell, consider for example, the finite linear 
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system of ordinary differential equations: 

^-u + A(t)u = f, u(0)=« , 



(1) 



where u is an unknown vector of ^components and A(t) is a given n x n matrix 
with bounded entries. In case A(t) is not continuous one does not expect u to 
be diffcrcntiable in the classical sense. Nonetheless, one would still like to make 
sense of the above initial value problem. To this end, one takes the Euclidean 
inner product, that we denote by (•, •), of equation (fTJ) with any vector field (j>(t), 
member of a family of conveniently smooth "test functions" that are identically 
equal to zero when t becomes large. Then formally (because no regularity on u 
is assumed) integrating by parts one can consider instead the equation 



A function u(t) is said to be a weak solution of (Q} if u satisfies equation 
© for every test function </> <G C°°([0, oo)), and 4>{t) is identically zero for t 
large enough. Notice that in order for u to satisfy ^ it does not need to be a 
differentiable function. 

In fact, without the refinement of functional analysis, this idea already goes 
back to Euler and Lagrange in the context of deriving the Euler-Lagrange equa- 
tion in the calculus of variations. Similarly, in the context of gas dynamics 
there is a good reason to introduce such a notion because compressible fluids 
support, and even when starting from smooth initial data may develop, shock 
waves which are discontinuous, and therefore have no derivative in the classical 
sense. Since Leray did not know much about the regularity of the solutions of 
the Navicr-Stokes equations it was crucial in his seminal work |57j to introduce 
weak solutions. Of course, this weak formulation leads to a loss of informa- 
tion which, for instance, may prevent the analysis of stability, or the proof of 
uniqueness. 

Weak solutions are defined as elements of some spaces of functions, which 
will be used below with no further introduction. 

With these prerequisites, a final and third thread will be the consideration 
of the equations of macroscopic perfect fluids and the relationship between them. 

We will start with the evolution of density, velocity and temperature (p, u, 0) 
given by the compressible Navier-Stokes-Fourier system of equations (cf. [S5] 
and [SD]): 



with the state equations for the pressure, p = p9, and the stress tensor: 




dtp + V • (pu) = , 

p(d t + u ■ V)u + Vp = eV • [zat(u)] 
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In these equations, for any vector u <G R d and any (scalar or vector valued 
function) v, the term (u ■ V)u denotes the expression 

i<^<d 

while ^ and k denote the kinematic viscosity and the thermal diffusivity of 
the gas, respectively. These quantities depend on the nature of the molecular 
interaction, while e denotes the dimensionless Knudsen number (proportional, 
up to a fudge constant, to the mean free path) which depends on the degree of 
rarefaction of the gas. A simplified system is then obtained (see, e.g., section 
1.1 in |60| ) for the fluctuations of order e of density, velocity and temperature 
about a constant background state, (p, 0,0), namely, 

p = ~p + ep, = + e6, u = eu. (3) 

In order to observe the cumulative effect of the small velocity, one has to consider 
the system for a time scale long enough (of the order of e _1 ), and hence by 
rescaling the time, by factor of e _1 , and following, for instance, the analysis 
presented at section 1.1 in [60], one obtains the following system: 



V-u = 0, V(0p + p0) = 0, 
~p(d t u + (u ■ V)u) + V7r = vAu , 

^p{d t + (u ■ V)0) = kA0, 

where it is the unknown pressure. This system is often referred to as the incom- 
pressible Navier-Stokes-Fourier system of equations. Observe that the equations 
for the velocity are independent of the temperature, while the advection term, 
(u- V)#, in the temperature equation depends on the velocity. The term incom- 
pressible comes from the V ■ u = condition, which implies that the Lagrangian 
flow map defined by: 

x i-> <f>(x, t) , <f>(x,0)=x, <fi(x, t) = u(4>(x, t), t) , 

is volume preserving. The equation V(0p + ~p0) = is known as the Bernoulli 
relation. Together with the equation for the temperature they describe how the 
flow, at this approximation level, propagate the initial fluctuations of density 
and temperature. In particular, if p = = the system reduces (hereafter the 
tilde will be omitted, p will be taken equal to 1, and following the tradition in 
this field we will replace the notation for the pressure tt by p) to the standard 
incompressible Navier-Stokes equations: 

d t u + (u ■ V)u - vAu + Vp = , V • u = . (4) 

In a typical engineering or physical situation characteristic length, velocity 
and time scales, L , U and T, respectively, are introduced. After rescaling the 
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relevant quantities in retaining the same notation for the rescaled dimension- 
less velocity, spatial and temporal variables, one obtains the undimensionlized 
Navier- Stokes equations 

d t u + (u ■ V)u - — Au + Vj? = , V ■ u = , (5) 
Re 

where 

Re=^ 
v 

is the dimensionless Reynolds number, which measures the ratio between a 
typical strength of the nonlinear advection term and a typical linear viscous 
dissipation in the Navicr-Stokes equations. Since equations ((H) and §5§ look the 
same (symbolically), we will abuse, hereafter, the notation and let v represent 
the inverse of the Reynolds number. 

In most practical applications the Reynolds number is very large (for in- 
stance, the Reynolds number is about 10 6 for the air around a moving bicycle, 
about 10 7 for a moving car, and is of the order of 10 12 in meteorological ap- 
plications). This naturally suggests to investigate and compare the behavior of 
the Navier-Stokes equations (O, for large values of the Reynolds numbers, to 
that of Euler equations, which are obtained formally by substituting Re = oo 
in ©: 

d t u + (u ■ V)u + Vp = , V • u = . (6) 

Thus, the third thread in our approach would be to compare the behavior of 
the solutions of (|SJ) to those of ©, as the Reynolds number Re — > oo. However, 
as it will be developed below, the relationship between these two systems is far 
from being simple. 

Since the above equations are derived for physical systems they are expected 
to respect certain global conservation laws which indeed can be deduced by 
simple computations, but only, assuming that the solutions are smooth enough 
to legitimize these computations. In particular, for sufficiently smooth solutions 
of the incompressible Navier-Stokes equations, say subject to periodic boundary 
conditions in a periodic box f2, by using the relation: 

/ ((u-V)u) -udx = f y2d x .(u j ^-)dx = 0, 
Jn Jn y z 

one obtains the formal energy balance: 
d f \u(x,t)\ 2 



dt 



dx + v \Vu{x,t)\ 2 dx = 0. (7) 



In the incompressible Euler and Navicr-Stokes equations the incomprcssibil- 
ity condition, V • u = 0, appears as a constraint on the dynamics of the fluid 
particles that are governed by the equation (Newton's second law): 

d 2 x(t) d 

= -u{x(t), t) = vAu(x(t), t) - Vp{x(t), t) , (8) 
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with the pressure, p, being the Lagrange multiplier corresponding to this con- 
straint. As such it can be eliminated by applying the Leray projection, P CT , 
to the equations (see, e.g., [55] for details concerning the Leray projection and 
other mathematical aspects of Navier-Stokes equations). For a domain without 
physical boundaries (i.e., the whole space or periodic box) P a is the orthogonal 
projection from L 2 onto the subspacc of divergence-free vector fields. That is, 
for any given w € L 2 , P a w is defined by solving the linear system: 

P a w = w + Vp and — Ap = V ■ w , (9) 

subject to period boundary conditions. In the presence of physical boundaries 
the Leray projection P a is the orthogonal projection from L 2 onto the subspace 
of divergence-free vector fields that are tangent to the boundary, hence P a is 
defined by solving system (|9|), subject to the Neumann boundary condition 

I = - <"» 

Eventually, due to the fact that our mathematical understanding of these 
problems and of their relations to turbulent flows is far from being complete, the 
use of explicit examples is important. Such an illuminating example is the shear 
flow solution of Euler equations, which is defined in fl = R 3 or in a periodic 
box fi = (R/Z) 3 by the formula: 

u(xx,X2,x 3 ,t) = (u 1 (x 2 ),0,u 3 (xi - tu 1 (x 2 ),x 2 ) ■ 

Indeed, it is a solenoidal vector field that solves Euler equations with the pressure 
p = 0. 

Other illustrative examples are the stationary solutions of Euler equations. 
Observe that, with q = p + Euler equations can also be written as 

d t u + (V Am) Au + \7q = , V • it = . 

Hence, stationary solutions are characterized by the relation (VAu)Am+V<7 = . 
In particular, any vector field that satisfies VAu = X(x)u is a stationary solution 
with q = (such solutions arc called Beltrami flows). In particular, any vector 
field that satisfies V A u = forms a stationary solution. In simply connected 
domains this relation implies that u is a potential flow, i.e., u = V</>, which is 
divergence-free, hence <p is a harmonic potential —A<p = . 

Finally, and for convenience we will denote by (•, •) the L 2 inner product in 
the relevant domain. 



3 The incompressible Euler equations 

In fact, as indicated by its name, this system of equations (one of the first 
systems of partial differential equations for fluid mechanics to be introduced) 
was derived by Euler in 1757, assuming the incomprcssibility of the fluid. It has a 
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very rich mathematical history with a very interesting recent important progress. 
The Euler equations in a domain ft (with or without physical boundaries) are 
given by 

d t u + (u ■ V)u + Vp = , V ■ u = , in ft. (11) 

If ft has physical boundaries then the system (|11[) is supplemented with the 
no-normal flow boundary condition: 

u ■ n = on 9ft . (12) 

Observe that the system (fTTj) corresponds to the limit case when the viscosity 
v = in or the Reynolds number is oo in ([5]). 

In the absence of viscosity (i.e., v = 0) any regular enough solution to the 
Euler system (fTTj) conserves the kinetic energy (compare with ([7]) for v = 0) . 

It turns out that the vorticity uj = V A u is "the basic quantity" both from 
the physical and the mathematical analysis points of view. The system (rj~T]) - (|12[) 
can be rewritten in an equivalent form in terms of the vorticity to give 

df io + u ■ Vcj = uj ■ V« , (13) 
V-u = 0,VAu = uj,u-ri = Oondfl. (14) 

The elliptic system given in equation (| 14|) fully determines u in term of uj (see, 
e.g., [35] and references therein), which makes the full system ([T3])-([TJ| "closed" . 
More precisely, the operator K : uj i-> u, defined by the unique solution to the 
elliptic system ([TT| . is a linear continuous map from C Q (ft) with values in 
C 1 '"^) (with a > 0), and from the Sobolev # s (ft) with values in # s+1 (ft) . 

Furthermore, for 2d flows the vorticity is perpendicular to the plane of mo- 
tion and therefore equation (jTB")) is reduced (this can also be checked directly) 
to the advection equation 

3 f w + !i-Vw = 0. (15) 

The quadratic nonlinearity in (|13[) has the following consequences, which we 
describe below, providing only the essential arguments, but omitting the details 
of the proofs. 

3.1 General analytical results for the 3d Euler equations 

The short time existence and uniqueness of classical solutions to the 3d incom- 
pressible Euler equations, with smooth enough initial data, was established (to 
the best of our knowledge) by Lichtenstein in 1925 [58] (see also [65], [66] and 
references therein for other proofs). 

Theorem 3.1 Let uq G C 1,Q (ft), with V • uq = and uq ■ n\gQ = 0. Then there 
exists a time T* > 0, that depends on uq, and functions u £ C ([— T*, T„], C 1,a (ft)) p| 
C 1 ([—T*,T*\, C°' a (Q)) , andp€ C°([-T t , T*], C*^ Q (ft)), such that the pair u,p 
solves Euler equations Hll\) - Uty) , with the initial data Uq. 
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In the proof of the above theorem, the "limit" on the lifespan of the solution 
is a result of the following nonlinear Gronwall estimate, that controls a certain 
norm, y(t), of the solution: 

y> < Cyi => y(t) < A . (i 6 ) 

Consequently, for a finite time, which depends on the size of the initial data y(0), 
the quantity y(t) is finite. To validate (|16|) and the above argument the initial 
data have to be chosen in the appropriate space of regular enough functions. In 
particular, if we consider the solution in the Sobolev space H s , with s > §, then 
by taking the scalar product, in H s , of the solution with the Eulcr equations 
(TTTj) . and using the relevant Sobolev estimates one obtains: 

Y^dh = -(("-^"."V* = -(v ■(«»«),«)*. (17) 

= ((«««) : Vu) Hs < C.ltiUir.lVtilio. < C\u\%. , 

which leads to (fT6|) and the local in time existence of a smooth solution. Observe 
that since V ■ u = we have used above the fact that (u ■ V)it = V • (it (g) it), 
where (it ® u) is the tensor (matrix) with entries a^- = itjUj, for i, j = 1, • • • , d. 

This notion of threshold of regularity for the initial data is well illustrated 
with the shear flow 

u(x,t) = {ui{x 2 ),0,u 3 (xi - tu 1 (x 2 ),x 2 ) , 

which is a well defined solution of Euler equations with pressure p = 0. Com- 
puting the derivative of 113 with respect to x 2 yields: 

td X2 ui(x 2 )d Xl U3{xi - tu 1 (x 2 ),x 2 ) ■ 

By investigating this term one observes that initial data which is in C 1,a generate 
solutions that will remain in C 1,a , for t > 0, and this is in full agreement with 
Theorem 13.11 of Lichtcnstcin. However, less regular initial data, say in C a , 
may, in some cases, generate solutions that will not remain in this space, for 
t > 0, but only in C°' a (cf. [9]). This counter-example shows that the Euler 
equations are ill-posed in C°' a , i.e. that Theorem 13. II of Lichtcnstcin cannot be 
extended to C°' a initial data. 

As it is often the situation in many nonlinear evolution problems, local reg- 
ularity means: local existence, uniqueness and local stability (i.e., continuous 
dependence on initial data). Moreover, one may exhibit a threshold for this ex- 
istence, uniqueness and propagation of regularity of the initial data (including 
analyticity Bardos and Benachour [5] , see also Kukavica and Vicol [55] and ref- 
erences therein). One of the simplest criteria for the regularity, up to the time 
T, is that the stress tensor Vu remains bounded in L 1 (0,T;L°°). To see that 
we observe that this bound implies the integrability of the Lipschitz constant 
of the vector field u, and hence the solutions of the ordinary differential equa- 
tion 4| = u(x, t) will never intersect and remains smooth within the interval 
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[0,T]. Consequently, and by virtue of ([TB"]) . the vorticity is transported, along 
the particles pathes, and it is stretched by Vu, and the latter is controlled by its 
L 1 (0,T;L°°) norm. However, to underline that it is the dynamics of vorticity 
which may induce singularities, a more sophisticated version of this criterion was 
established by Bcale, Kato and Majda [T2], to only retain the anti-symmetric 
part of the stress tensor (Vu — (Vu) T ) (equivalent to the vorticity vector ui); 

Theorem 3.2 Beale-Kato-Majda flf Let u(t) e C([0, T),H s )f] (^([O, T), 
for s > | , be a regular solution of the 3d incompressible Euler equations in the 
interval [0, T). Assume that 

[ |VAu(.,t)|L«di<oo, (18) 
Jo 

then u(t) can be uniquely extended to the interval [0, T + S], for some S > 0, as 
a regular solution of the Euler equations. 

While a companion result of Constantin, Fefferman and Majda [55] to The- 
orem [3T2] states that it is mostly the variations in direction of the vorticity that 
may produce singularities. 

Theorem 3.3 [25] Let u be a smooth solution of the 3d Euler equations, which 
is defined in Q = Q X [0,T). Introduce the following quantities (which are well 
defined for every t £ [0, T) ): 

fci(i) = sup|u(x,t)|, (19) 

which measures the size of the velocity, and 

k 2 (t)= sup 1 : , (20) 

x,yen,x^ v \x-y\ 

which measures the Lipschitz regularity of the direction 

e(*.*) = i^ 

\uj(x,t)\ 

of the vorticity vector. Then under the hypotheses: 

(fei(t) + k 2 {t))dt < oo and / k 1 (t)k 2 (t)dt < oo (21) 



the solution of the 3d Euler equations exists, and is as smooth as the initial 
data, on the interval [0, T + S], for some S > . 

It is worth mentioning that a similar criterion to Theorem l3.2[ involving the 
symmetric part of the stress tensor S(u) = i(Vu+(Vu) T ), was introduced later 
by Ponce [72] • See also [34] for other related blow-up criteria for the 3d Euler 
equations. Moreover, see [46] and references therein for a recent survey of the 
blow-up problem. 



3.1.1 About the two-dimensional Euler equations 

In the 2d case, the vorticity lu = V A u obeys the evolution equation 

dtUJ + u- Vcj = 0. (22) 

This transport evolution equation enforces the conservation of any LP norm 
(1 < P < oo ) of the vorticity vector. Taking advantage of this observation in 
1963 Youdovitch proved in his remarkable paper [5U] the global existence and 
uniqueness for all solutions with initial vorticity in L°° . If the vorticity is only 
in LP for 1 < p < oo one can prove the global existence (without uniqueness) of 
weak solutions. The same result holds also for p = 1, and also for the case when 
the vorticity is a finite Radon measure with definite sign (cf. Delort [33]) or 
"simple" changes of sign (cf. [62] and [64]; see also [7]). It is worth mentioning 
that the proof is more delicate in these limit cases when p = 1 or when w is a 
bounded Radon measure. 



3.2 Weak solutions of Euler equations 

Definition 3.1 A weak solution, on the time interval [0,T], of the Cauchy 
problem for the Euler equations ill]) , with divergence-free initial data u(x,0) £E 
(L 2 (EL d )) d , is a divergence-free vector field 

u(x, t) € L 2 (0, T; R d ) d p| C([0, T}; (L 2 (R d )) d w $ 

and a distribution p which, for any smooth test function (j>(x, i) vanishing for t 
near T , satisfy the Euler equations 111]) (or ^) in the distribution sense, that 
is: 

\ [ u(x,t) ■ d t (j){x, t)+(u(x, t) <g) u(x, t)) : ^ x 4>{x, t) ] dxdt+ (p, V x 4>) 

JR d 

u(x, 0) • 4>{x, 0) dx , 



(23) 

where (•,•) denotes the distributional action. 

It is important to observe that the hypothesis u(x, t) £ L 2 (0, T; R d ) d is sufficient 
to insure that the matrix u ® u is well defined as a distribution (i.e., no need 
for any additional hypotheses on the partial derivatives of u). 
Moreover, equality (|23j) is equivalent to the following: 

/ [ u(x,t) ■ d t <j){x,t) + {P IJ (u(x,t) <g) u{x,t)) : V x <f>(x,t)) } dxdt 

(24) 

u{x, 0) • 4>(x, 0)dx , 



1 The symbol w refers to the weak continuity. In the present setting this means that for 
any (f> € L 2 (R d ) the functiont i-> (u(t) , phi(t)) is continuous. Moreover such property follows 
from the relation 1231 1 
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with P a denoting the Leray projection given in ©. Eventually, this is also 
equivalent to assuming that the relation 




I / [ u(x, t) ■ dt<j){x, t))+(u(x, t) ® u(x, t)) : V x <j>{x, t) ] dxdt 
o Jn d 



(25) 




u(x, 0) • <p(x, 0) dx 



holds for every divergence-free test functions cj>. 

Moreover, there are other variants of the notion of weak solutions. One 
of the most efficient, introduced by DiPerna and Lions [50], is the notion of 
dissipative solution. It is defined in terms of the stability with respect to any 
smooth divergence-free vector field w (which may not be a solution of the Euler 
equations) . Such a smooth vector field w evidently satisfies the relation: 



where E(w) is the residual, satisfying P a E = E. Notice that E = if and only 
if iy is a smooth solution of the Euler equations ©. Next, observe that if u is 
a smooth divergence-free vector field then by integration by parts, with respect 
to the spatial variable, one obtains: 



If in addition we assume that u is a smooth solution of the Euler equations then 
from ([26]). ([6]) and the relevant boundary conditions (e.g., periodic boundary 
conditions or (fT2]) ) we conclude : 



dtW + P a {{w V)w) =E, 



(26) 




with 





which, upon integration with respect to time, gives 



2 / / \E(w(x, s)) ■ (u(x, s) — w(x, s))| dx ds 




+ 2 / / \(u(x, s) — w(x, s))S(w)(u(x, s) — w(x, s))\ dx ds 
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and by virtue of Gronwall inequality yields 

\u(x,t) -w(x,t)\ 2 dx < ( [ \u(x,0)- w(x,0)\ 2 dx)e 2 ti l s ( TO )H-W 
^ Jn ' 

2 [ ( j |S(w(a;,s))-( W ( : E,s)-u;(x,s))|da;)e 2 ^ l|s(,u)ll ~ (r)dT ds 



(27) 



Inspired by the above we introduce the following definition of dissipative solution 
for Euler equations (j5]): 

Definition 3.2 A dissipative solution of Euler equations (0|) is a divergence- 
free vector field u which together with any smooth divergence- free vector field w 
(which need not be a solution of Euler equations) they satisfy the inequality: 



n 

t 



\u(x,t)-w(x,t)\ 2 dx < U \u(x,0)-w{x,Q)\ 2 dx S j e a/o*|S(«)llo.««to 

n (28) 
2 / \(d t w(s) + P a ((w(s) ■ V)w(s)),u(s) - w(s))\ e 2 ^ H^MIUM^ ds 



n 



Remark 3.1 With the above definition one should keep in mind the following 
facts. 

1. Inequality V28\) makes sense with no other regularity conditions on u other 
than the fact that it belongs to (L 2 ([0,T; (L 2 {Vi)) d )) d (d = 2 or d = 3). 

2. If w is a genuine smooth solution to Euler equations then one has E(w) 
and inequality i28\) reduces to 



\u(x,t) ~ w{x,t)\ 2 dx < (f \u{x,Q)-w{x,Q)\ 2 dx^ e 2 J? 1^)11 -(«)««« , (29) 

which implies both uniqueness of dissipative solutions, within the class of smooth 
solutions, and stability of dissipative solutions with respect to smooth solutions. 

Below it will be observed that weak solutions of Euler equations © may 
not conserve energy. Therefore, the connection between weak and dissipative 
solutions of ((6|) requires some additional information on the control of energy. 
Following [29] we introduce the following: 

Definition 3.3 A divergence-free vector field u £ L 2 (0,T;R d ) d is called an 
admissible weak solution if it is a weak solution and if it satisfies the additional 
hypothesis: 

/ \u(x,t)\ 2 dx < / \u(x,0)\ 2 dx, for all t € [0, T] . (30) 
Jn Jn 

As a result of the above definitions we have: 

Proposition 3.1 Let u be an admissible weak solution of Euler equations (0), 
in £1 = H d or in a periodic domain ft = (R/Z) d ; then u is a dissipative solution. 
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Proof Below, we only present the main idea of the proof. Let w be any smooth 
divergence-free vector field. The point is to compare 

— f \u(x, t)— w(x, t)\ 2 dx — — J \u(x, t)\ 2 dx+ — J \w(x,t)\ 2 dx— f (u(x,t)-w(x,t))dx 
2 Jo 2 J n 2 Jn J n 

and 

- [ \u{x,0yw{x,0)\ 2 dx = i / \u{x,Q)\ 2 dx+l- [ \w(x, 0)\ 2 dx- j (u(x,0)-w(x,0))dx 
2 2 Jq 2 Jq Jq 

With hypothesis (f3T)l ) it is then enough to show that: 

i/" \w(x,t)\ 2 - J (u(x,t) ■ w(x,t))dx = i /" |w(x,0)| 2 - y (u(x,0) • w(a;,0))dx 
+ y — ^— y s)| 2 — y s) ■ w(x, s)) dx^J ds 



< - I \w(x, 0)| z - / (u(x, 0) ■ w(x, 0)) 



o Jo 
t 



(u — — w) dxds+ 

E{w(x 1 s)) ■ (u(x, s) — w{x, s)) dxds , 

/0 JO 

or that 

E(w(x, s)) ■ (u(x, s) — w(x, s)) dx . 

Here, the most nontrivial step that needs further justification is the one involving 
the term J„ (d s u(x, s) ■ w(x, s)) dx. But since w is smooth vector field the notion 
of weak solution can be invoked. 

However, for Eulcr equations the introduction of weak solutions corresponds 
to an important loss of information illustrated by the following theorems of De 
Lcllis, Szekelyhidi and Wiedemann. 

Theorem 3.4 1. Let f2 C R d , for d = 2,3, be any bounded open set, and 
let T > 0. Then for any e(x,t) > 0, continuous function defined on Q = 
D, x [0, T] C R d x Rt, there exists an infinite set of weak solutions of Euler 
equations in R rf which satisfy: 



\u{x,t)\ 2 



e(x,t), for all(x,t) <S Q, and \u(x,t)\ = 0, for all(x,t) ^ Q. 



2. There is an infinite set of initial data uq € L 2 (R rf ) for which there exist 
admissible solutions (hence they are dissipative) . 

3. For any initial data uq € L 2 (R rf ) there exists infinitely many weak solutions 
u(x,t) £ C(Rt; L 2 (R, d ) w ), such that u(x,0) = «o- 
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The first two statements in Theorem l3.4l are due to De Lellis and Szekelyhidi, 
and the last one is due to Wiedemann following their approach. A complete 
proof of the above theorem is out of the scope of this contribution, but in the 
following subsection we give some comments which might be appropriate. 

Comments on Theorem 13.41 The first two statements of Theorem 13.41 
have, in fact, a long history over the second half of the previous century. In 
1993 Schcffer [76] was first person to observe the existence of weak solutions of 
Euler equations with compact support in time. Later, in 1997, Shnirelman [ST] 
refined Scheffer's proof. The basic method was to construct approximate highly 
oscillating (in space and time) divergence-free vector fields Uk (with compact 
support in space and time), that would converge strongly in L 2 (R d x R t ) (in 
order to be able to pass to the limit in the nonlinear term Uk £3> Uk), while by 
equation ([26]) the sequence of residuals: 

E k := d t u k + P CT ((V ■ K ® u k )) , 

would converge weakly to zero. 

Of course, such weak solutions, for which De Lellis and Szekelyhidi coined 
the name "wild", are not physical at all. They describe a fluid flow which 
starting from rest and without any external force applied to it, begins to move 
and then goes back to rest. If this situation could be applied to "real world fluid 
flows " that would solve the energy crisis! 

To establish their proof, De Lellis and Szekelyhidi realized a striking anal- 
ogy between "wild solutions" of Euler equations and the classical problem of 
isometric imbedding. 

The simplest example of such an imbedding problem goes as follows: 

Let S 1 — {x G R 2 : |a;| = 1} C R 2 be the unit circle, and, for any r > 
0, B 2 (r) — {x G R 2 : |x| < r} C R 2 be the closed disc of radius r. The 
question is: what is the best regularity (in term of the Holder exponent /3) 
that is independent of r, and for which there exists a map <f) r <G C^(S 1 ; B 2 (r)) 
(dubbed isometric imbedding) which enjoys the following property: <f> r {S 1 ) is, 
for the Riemannian metric induced by R 2 , a manifold isometric to S 1 ! 

Now if one requires <fi r G C 2 , the curvature of <f> r (S l ) should be bounded 
from above and below, independent of r. While to place the curve ^(S 1 ) in a 
ball of arbitrarily small radius, one should expect its curvature to go to infinity, 
asr->0. This is an obstruction which was realized by Gauss under the name 
of "egregious theorem". Hence, to realize the isometric imbedding, one has 
to relax the regularity of the mapping cf> r from being C 2 . What was proven 
(cf. [30] and the reference therein) is the existence of a small enough Holder 
exponent a G (0, 1), which is r independent, such that for any r > there 
exists an isometric imbedding <j> r G C l a of S 1 into ^(S 1 ) C B 2 (r). The main 
idea of the proof is to show first the existence of a mapping tp r that reduces 
the distance, on the Riemannian manifold ^(S 1 ), and then accumulate (add) 
small oscillations (wiggles of the curve) to reach the equality of the lengths; and 
thus prove that the isometric imbedding can be realized with the limit function 
of this process, 4> r G C 1,a . 



14 



In the same way one can rewrite Euler equations as a linear system with a 
nonlinear constraint, namely, it is equivalent to the system: 



d t u + V x -M + Vq = 0, u(x,0) = u (x) ,V • u = 0, (31) 
with the constraint 

\u\ 2 

M = u®u- -j-I, (32) 

where M is a traceless self-adjoint matrix, and / is the identity matrix. To solve 
the above system one considers, instead, the following more relaxed problem. 
Let Q and e be as in Thcorem l3.4l Let A be is the closed hull (in some convenient 
topology) of the set of all triples (it, q, M), where u is a vector field, q is a scalar, 
and M is a self-adjoint traceless matrix, which satisfy system (|3"Tj) . and the 
following relaxed constraint: 

2e(x t) 

(u®u)-M < — -r^I, for all (x, t) e Q, and \u(x, t)\ = , for all (x, t) $ Q. 



The elements of A are called sub-solutions of (|31l) - (|32|) . Then the proof of 
Theorem 13.41 goes as follows. First, one should show that the set of sub- 
solutions, A, is not empty. Then, to conclude the proof, one constructs a 
sequence (uk,q k , Mk) £ A, which is constructed by adding, more and more, 
spatial and temporal oscillations of small amplitudes and high frequencies, and 
which converges strongly, in L 2 , such that: 



lim 

k— >oo 



(u k ®u k )-M k -^I 



. (33) 



Remark 3.2 1. Uniqueness among dissipative solutions, whenever one of them 
is regular, implies that the wild solutions stated in point 3 of Theorem \3.4\ (for 
every initial data uq € L 2 ) cannot be admissible whenever the initial data is 
regular. Indeed, as Wiedemann observes in |#7| / his construction of sub-solutions 
requires a jump in energy at time + . 

2. Similarly the initial data giving rise in point 2 of Theorem \3.4\ to admissible 
solutions cannot be regular; otherwise, the solutions being admissible would co- 
incide with the unique dissipative solution. On the other hand, since the proof 
involves the Baire category theorem, the construction of these initial data is, in 
general, far from being explicit. However, for the case of Kelvin- Helmholtz type 
initial data: 

u (xi,x 2 ) = {-sign(x 2 ),0) , (34) 

in a periodic domain, explicit sub-solutions, and then infinitely many solu- 
tions, have been provided in i82f . In 3d, and for initial data of the form 
(—sign(x2),Q,w(xi,X3), the construction presented in J82f gives rise to an infi- 
nite set of admissible solutions, cf. fWjj. 

3. Eventually, with a better adapted choice of the high-frequency corrections 
to the sub-solutions, one can improve the spatial and temporal regularity of the 
wild solutions obtained at the limit. Along these lines, using Beltrami flows as 
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building blocks for the high-frequency corrections in one proves, for every 

fixed a in (0,1/10), the existence of admissible wild solutions which belong to 
C 0,a , and which exhibit strict energy decay. 



3.3 Energy decay and turbulence 

For then Euler equations the only natural estimate is the energy estimate which, 
as mentioned above, gives for smooth solution the relation: 

f \u(x,t)\ 2 dx = - [ \u{x,0)\ 2 dx. (35) 



2 j , v > 2 

From this relation one deduces that any e dependent family u e of smooth solu- 
tions is uniformly bounded according to the formula: 

Vt, e \u e (t)\ L 2 < sup \u e (0)\ L 2 < oo . (36) 

e 

Variants of this energy estimate (with inequality instead of equality) can be 
obtained and are valid in most of the physical situations (including Navier- 
Stokes with weakly converging initial data and viscosity going to together 
with boundary effects) . 

This implies the existence of subsequences, still denoted by u e , which are 
weakly converging to a limit (lim e _>.o u e )(x,t) e i°°(R f ;L 2 ) . However, at this 
level of generality it is not possible to conclude that this limit remains a solution 
of Euler equations. The problem is that the weak limit of the tensor product 
may not be the product of the weak limits: 

lim(u e ® u e ) ^ (lim u e ) ® (lim u e ) (37) 



and therefore the weak limit u = lim e _j.o u c is solution of the equation: 

(38) 



V-u = 0, 9 t u + V(u(g)M) + i?T + Vp = 0, 
with RT = lim(M e — u) <E> lim(u £ — u) . 



One simple example (which is due to DiPerna and Majda [35]) can be con- 
structed with the shear flow: 

uAx, t) = (sin( — ), 0, w(xi — t sin( — )) 
e e 

which is uniformly bounded and which weakly converges to 

It = (0,0, / w(xi — t sin s)ds) . 
Jo 

Obviously this is not a solution of the Euler equations, by direct inspection and 
also because the term u\d Xl w e does not go to zero! 

Weak limits can be interpreted in term of "coarse graining", as the deter- 
ministic counter-part of the Reynolds stress tensor for small fluctuations in the 
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statistical theory of turbulence. In this setting one could consider, as turbulent 
flows which are not solutions of the Euler equations due to the appearance of 
the extra term RT in (|38[) . Now the formula 

for almost every t lim \u(t) - u e (t)\ 2 = lim \u e (t)\ 2 - \u(t)\ 2 (39) 

e— >0 e— >0 

implies that the conservation of energy would give RT = and therefore no 
turbulence. Hence, this implies the conjecture that turbulence would be related 
to anomalous energy dissipation. 

Indeed, the issue of the conservation of energy for weak solutions was already 
raised by Onsager |71) when he formally wrote: 

0= / " [V • (u ® u)] ■ udx ~ y*(V3 U • V^u • V^u) ~ J\\7^u\ 3 dx. (40) 

He conjectured that any weak solution that would be more regular than C 0, 3 
would preserve energy. A rigorous proof of this result, inspired by the naive 
formula (|40|) . was given by Constantin, E and Titi [24]. Of course this does 
not mean that the conservation of energy implies the regularity of the solution; 
simple examples of singular weak solutions which preserve energy can be con- 
structed with the shear flow [9] or, as noticed above following ( 30 ), as "wild" 
solutions. Eventually to complete this picture it would be interesting to exhibit 
weak solutions less regular than C°'3 for which energy would decay 

The first step was the construction of continuous (in space and time) solu- 
tions with decaying energy (cf. [3 2) for the spatially three-dimensional case, and 
[2"0] for the spatially two-dimensional case.) One of the basic ingredients in the 
spatially three-dimensional case was the introduction of the Beltrami flows in 
the construction. Then a refinement of the above method gave the existence of 
energy decaying weak solutions which belong to C°' Q , with a < ^ (cf. [31] for 
the spatially three-dimensioanl case, and |19j for the spatially two-dimensional 
case). Most recent improvement, due to Isett [50], extend the Holder regularity 
of these energy decaying solutions to a < in a preprint of 171 pages! 

4 Navier- Stokes equations without and with bound- 
ary effects 

4.1 Why in some cases Navier-Stokes rather than Euler 

In the section[2]thc Euler equations were presented as a formal limit, as Re — > oo, 
of Navier-Stokes equations. However, we have also added the the relation be- 
tween these two equations "is far from being simple." The Euler equations were 
derived first, and it took about 45 years to obtain Navier-Stokes equations as 
an "improved" model (cf. Navier [5§] and Stokes, [ED]), taking into account the 
viscosity effects which are better adapted to the real fluids . Of course, at that 
time this was not motivated by the arguments of De Lellis and Skzekelyhidi 
on wild solutions, but by the need for a mathematical description of the fluid's 
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viscosity. Moreover, the introduction of the viscous term —vAu, with the corre- 
sponding boundary conditions, is essential to describe the interaction of a fluid 
flow with physical obstacles. In general, this physically observed interaction 
cannot be described only by means of Euler equations, subject to the imper- 
meability (no-normal flow) boundary conditions, as it will be illustrated by the 
following proposition: 

Proposition 4.1 Consider an obstacle K C R 3 , with a smooth boundary which 
is topologically isomorphic to a closed ball. Let u be the velocity field of a sta- 
tionary solution of Euler equations in fl — R 3 \if , that is constant at infinity, 
i.e., limui^oo u(x) = u^, and which is satisfying the impermeability (no-normal 
flow) boundary condition, u ■ ft — on dQ. Then the total force exerted by the 
fluid flow on the obstacle K is 0. 

Proof (cf. [H] or for details) Since the fluid flow is stationary and 
constant at infinity it is also irrotational and therefore one has the Bernoulli 
theorem: 

V^+p)=0. (41) 

Moreover, since the obstacle K is homcomorphic to a ball in R 3 the domain 
f2 = H 3 \K is simply connected; therefore the velocity field u fluid is both 
"potential" and solenoidal . Hence, it can be described by the system: 

u = -V<t>, -A0 = O (42) 

Thus, with the asymptotic behavior, at infinity, of the Green function of the 
Laplacian in R 3 (behaving like (47r) _1 |2;| _1 ) one has: 

u{x) = Woo + 0(M~ 3 ) and p(x) = Poo + 0(\x\~ 3 ) , (43) 

where the above asymptotic behavior of the pressure follows from (|41l) . Now, 
to compute the force exerted by the fluid one introduces the ball Br = {x £ 
R 3 , |a;| < R}. For R large enough one has K C Br. Therefore, thanks to the 
impermeability boundary condition, u ■ ft = 0, one can compute, F, the force 
exerted on K, as follow: 



pftda = — I (it ■ ft)u + pft da 
on Jan \ J 

/ («• f )u+P^ ]da- / ((«•■! 

J\x\=R \ K K ) J\x\=R \ K 

/ f (u ■ V)u + Wp)dx- ( 

JnnB n \ J J\x\=R \ 

J\x\=R \ U K J 



X X \ 



(44) 
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where the last equality is due to (|4T|) . Finally, by virtue of (|4U| the relation (|44p 
implies: 

F = - I piida = - lim / [ (u ■ — )u + p— I da = ! (45) 
Van y|a;|=fl V R R ) 



Remark 4.1 The above computation was already known by d'Alembert, for this 
reason it is called the d'Alembert paradox. A similar computation can be per- 
formed in two dimensions, but in this case the domain Q = Tt?\K , where K 
is isomorphic to a disc in R 2 is not simply connected. Therefore, irrotational 
vector fields need not to be potential. In addition, the Green function in two 
dimensions behaves like ^logj^j, and it produces a nonzero contribution (lift) 
called the Kutta-Joukowski relation. 

In conclusion, we observe that the introduction of the viscosity, and of the 
corresponding boundary condition associated with to the Laplacian, may lead 
to effects that cannot be described by Euler equations (which corresponds to 
the singular limit v = 0), and that these effects might become important for 
"very small" values of v , where v > 0. 

For example, the 2d Navier-Stokcs equations in the half-plane, X2 > 0, sub- 
ject to Dirichlet boundary condition, is a well-posed problem; and its vorticity, 
u), evolves according to the equation 

dtOJ v — vAujv + u v ■ Vw„ = . (46) 

In the inviscid case, v = 0, the vorticity, ui, evolves according to: 

d t oj + (it • V)w = , (47) 

which does not indicate any interaction of the vorticity inside the domain with 
the boundary R x {x2 = 0}. However, for the viscous case, v > 0, the values of 
the vorticity at the boundary, R x {x2 = 0}, will depend in a nontrivial way on 
the pressure . To see that, we compute the normal derivative of the vorticity 
at R x {x2 = 0}, and use the incompressibility condition, V • u v = 0, and the 
Dirichlet boundary condition, u v = at R x {x2 = 0}, to obtain: 

d X2 u v = d X2 (d Xl ul-d X2 ul) = d Xl (d X2 ul)~dl 2 ul = -d^ul-d^ul = -d 2 X2 u\ . 

Next, use the evolution equation for the first component of the velocity, u x v , in 
the Navier-Stokes equations at the boundary, R x \x2 =0}, to get: 

vdl 2 ui = d Xl p. 

Combining the above yields the so called Lighthill equation: 

- vd X2 u v = d Xl p, atR x {x 2 = 0} . (48) 
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The above relation, (|4"g)) . implies that the variation of pressure at the boundary 
generates the vorticity there or vice versa. Notice that the quantities involved 
in (|48p may be very large, for small values of v > 0, leading (cf. section (|4.4[) ^ 
to important effects and interaction with the boundary. 

To measure such effects, in particular in comparison with the Reynolds num- 
ber, it is convenient to consider a slightly more general class of boundary con- 
ditions. Therefore, in an open smooth domain O C R d , with d = 2 or d = 3, 
one supplements the Navier-Stokes equations 

d t u u + {u v ■ V)u v — v£±u v + Y7p„ = , in ft (49) 

with the boundary conditions: 

u v ■ n — , v{dnU v + C(x)u v ) T + Ait„ = on dfl , (50) 
with A(i/, x) > , C{x) e C(R"; R n ) , (51) 

where n is denoting the outward normal to <9f2, and a T = a — (a ■ n)n is the 
tangential component of the vector a at the boundary d£l. Of course, when the 
problem is considered in the whole space, or in a periodic box, these boundary 
conditions are omitted. In (|50p there are two statements which represent, on one 
hand, the impermeability condition u v ■ ft = , which is in some sense universal 
(independent of the fluid, and in particular of whether it is viscous or inviscid) 
and, on the other hand, a condition which models the friction between the fluid 
flow and the boundary, and which may exhibit some nontrivial behavior in terms 
of the Reynolds number. With \(v) = oo, the second statement corresponds to 
the Dirichlet boundary condition: the component of the velocity on the plane 
tangent to dfl being also zero and, therefore, it is called the "no-slip" boundary 
condition. 

In the more general case, C(x) is a linear bounded operator acting on the 
space of tangent vectors to dSl with the L 2 (d£l)— norm. With a convenient 
choice of C(x) this statement includes the following conditions 

(S(u„) ■ n) T = {dau v ) T - (V*n • u v ) T => v{S{u v ) ■ n) T + \u v = , (52) 

and 

(V A «„) A n = {dnu v ) T + (V*n • u„) T ^> u{V A u v ) An + \u v = . (53) 

Condition (|52[) is the standard Fourier law, which can be either derived by 
phcnomcnological arguments, or deduced from a boundary condition for the 
Boltzmann equations, due to the atomic structure of the wall (cf. pQ, [6] and 
[68]). Condition ([53| emphasizes the role of the vorticity. 

Integrating by parts after multiplying equation (|49[) by u v , and using some 
standard trace estimates one obtains the: 

Proposition 4.2 Any smooth solution of the Navier-Stokes equations, 
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with the boundary conditions i5U\) , satisfies the relation (energy balance): 



4.2 Basic mathematical results for the Navier- Stokes equa- 



From the estimate (|54|) one deduces, and this was the original contribution of 
Leray [57], that for any initial data uq € L 2 there exists (globally in time) 
weak solution to the Cauchy problem. However, there are several important 
differences between the 2d and 3d cases. In the 2d case one can prove that for 
any uq € L 2 the solutions are unique and that the formal energy equality (|54p 
is exact and can be justified rigorously. In the 3c? case the situation (which is an 
object of the present Clay prize) is more complex. With no restriction on the 
size of the viscosity, v, when it is compared inversely to the spatial oscillation 
in the initial data (even very smooth initial data), the only result available 
is the (global in time) existence of a weak solution which satisfies the energy 
inequality, i.e. instead of the equality sign = in (|54l) one has the inequality 
sign <. Such weak solutions that satisfy the energy inequality in ((54)) are called 
Leray-Hopf solutions. It is worth mentioning that all the known existence proofs 
of weak solutions to the Navier-Stokes equations lead at the end to Leray-Hopf 
solutions. The uniqueness of weak solutions, or Leray-Hopf weak solutions, and 
the propagation of regularity, for any time, of the initial data are still open 
problems. 

There is a long history of question of global regularity and of blow up criteria 
of smooth solutions of the 3d Navier-Stokes equations, taking a full advantage 
of the regularizing effect of the Laplacian. As observed in Theorem 13.21 the 
absence of blow up for the solutions of the 3c? Euler equations is enforced by 
the Bealc-Kato-Majda condition, which concerns the vorticity. For the Navier- 
Stokes equations it is enough to assume that the velocity remains bounded to 
ensure the regularity. This condition can be widely improved starting from the 
pioneer works of Prodi [74] and Serrin [78], involving the L p ((0, T), L 9 (f2)) norms 
of the velocity field, with ^ + ^ = 1 (for more details see, for example, the survey 

2 The appearance of unbounded, or very large velocity, is not in inconsistent with the 
physical derivation of the equations. This is because, as it has been mentioned in the "Pre- 
liminaries" section (section O , the incompressible Navier-Stoke equations can be deduced as 
a small Mach number approximation. In this case the velocity u e is replaced by eS e , and the 
fact that u t is bounded does not imply that u e , the solution of the "incompressible" equations, 
needs to be also uniformly bounded. 




(54) 



and the v— uniform estimate: 




(55) 



tions 
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articles [35] [53], [53] and references therein). Most recently, there has also been 
some criteria for global regularity of the 3d Navier-Stokes equations involving 
the pressure (cf. [T7] and references therein). In addition, some other sufficient 
regularity conditions were recently established in terms of only one component 
of the stress tensor Vw, (cf. [TB] and references therein). Furthermore, in the 
presence of certain spatial symmetries, such as axi-symmetric [55J or helical [53] 
flows, the Navier-Stokes equations are known to be globally well-posed, but such 
results are still out of reach for Euler equations (see, e.g., [37] and references 
therein for some related results concerning helical Euler equations). It is worth 
mentioning that vorticity stretching term, which is the main characteristic of 
'id flows verses 2d flows, is nontrivial for axi-symmetric and helical flows. 

Eventually, the two following sharp results (even, if once someone proves that 
weak solutions of Navier-Stokes equations are always smooth) have a genuine 
interest by the detail and the physical interpretation of their proofs. 

• Caffarclli, Kohn and Nirenberg |16j : 

For any initial data uq € L 2 there is a weak solution (called a suitable 
solution) for which the support of its singularity set, in space-time, is contained 
in a set whose one-dimensional Hausdorff measure is zero. This means that the 
singularity set, in space-time, is "less" than one-dimensional. 

• Seregin and Sverak [77] : 

Smooth solutions of the 3d Navier-Stokes equations remain smooth for as 

long as either the pressure is bounded from below, or the Bernoulli pressure 
i i 2 

— h p is bounded from above. 
For instance what [T5J says is that the velocity of the fluid cannot be too big, 
say |u| > M, except in regions which are small with respect to M. The result of 
[77] can also be interpreted by saying that for as long as there are no cavities in 
the flow the solution can be uniformly controlled. Of course, the minute cavities 
are formed the Navier-Stokes equations are no longer valid as a physical model. 

4.3 Zero viscosity limit of solutions of Navier-Stokes in 
the absence of boundary 

As mentioned in the previous section, the introduction of the viscous term, i.e., 
the Navier-Stokes equations, is essential in order to take into account boundary 
effects. Conversely, in the absence of physical boundaries, one has the following 
proposition: the behavior, as the viscosity v — > 0, of Leray-Hopf weak solutions 
of the Navier-Stokes equations is related to the behavior of the solutions of the 
Euler equations according to 

Proposition 4.3 1. For given initial data, u 6 L 2 (R d ), Leray-Hopf weak 
solutions of the Navier-Stokes equations are bounded ( uniformly with respect to 
v) in L°°(R. t ;L 2 {Ri)) . 

2. Any weak limit, u, as v — > 0, of such a family of Leray-Hopf weak solutions 
of Navier-Stokes equations is a dissipative solution of Euler equations. 

3. Assume that, with uq as initial data, Euler equations have a smooth solution 
u{t), on the time interval [0, T]. Then u(t) = u(t), for all t G [0,T]. 
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Remark 4.2 The issue in the above proposition is that, in the absence of bound- 
ary effects (with boundary the situation is completely different, see section \4^4\ 
below), any weak limit of solutions to the Navier-Stokes equations is a dissi- 
pative solution to the Euler equation. Then the next issue does not concern 
Navier-Stokes equations but the dissipative solution of the Euler equations. In 
the presence of a smooth solution (with the same initial data) of Euler equations 
point 2 of Poposition \3.1\ implies that this weak limit is the smooth solution of 
the Euler equations. Hence, it is unique and nothing spectacular happens as 
described with more details in the RemarkU.3[ below. 



Proof of the proposition 14.31 : By definition, a Leray-Hopf weak solutionsb 
u v {x,t) to Navier-Stokes equations satisfies the energy inequality (here without 
boundary terms, since there is no physical boundary): 



\ ( \u u (x,t)\ 2 dx + p [ [ \Vu„(x,s)\ 2 dxds < ^ [ \u a (x)\ 2 dx . (56) 

2 Jn d Jo JTL d 2 JTL d 

Therefore, u v is uniformly bounded in L°°(R t ; L 2 (R^)) and, modulo extraction 
of a subsequence, converges weak—* in this space to a limit u. Proceeding as 
in section 13.21 (for rigorous justification of all the steps below one can follow 
the arguments introduced in |79|). for any divergence- free smooth test function 
w(x, t) satisfying 

dtw + P<r((w-V)w)=E(w), (57) 



one obtains the estimate: 



\u v {x, t) — w(x, t)\ 2 dx + 2u 



R'< 



/ / |Vu„(x, s)\ 2 dxds < 
Jo Jn d 



2 / / \(E(w(x,s)) • (u(x,s) — w(x,s)))\dxds +2u / / (Vii„(a;, s) • Vw(i, s))dxds 
Jo iR i Jo Jn d 

+ 2 / / \(u v (x, s) — w(x, s)S(w)(u v (x, s) — w(x, s))\dxds 

Jo JR d 



\uq(x) — w(x, 0)\ 2 dx . 



Thanks to Grdnwall Lemma one has: 



(58) 



R'' 



+ 2 



+ 2v 



\u v (x,t)-w(x,t)\ 2 < f / \u (x)-w(x,0)\ 2 dx le 2 ^ |S(«)|U(«)<ta 



R'' 



R'< 



\(E(w(x,s)) ■ {u„(x,s)-w(x,s)))\e 2 ^ ls{w)l ^ T)dT dx 



ds (59) 



R-' 



(yUu(x, s) ■ Vw(x, s))dxe 2 SI |S(«-)Hoo(r)dr 



ds . 
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Thus, by virtue of the energy estimate: 

2v I I \Vuv(x, s)\ 2 dxds < / \uo(x)\ 2 dx 
Jo Jn d Jn d 

there exists a subsequence of u v , which will also be labeled u u , such that ^JvSIu v 
converges weakly in L 2 ([0, T]; L 2 (R^)). As a result one has: 



lim 2v 

v-yO 



R'' 



(V« w (x, s) ■ Vw(x, s))dx e 2 SI |S(«0IL(r)dr 



ds = . 



Moreover, as noticed above, the weak convergence gives 



\u(x, t) — w(x, t)\ 2 < liminf / \u v (x, t) — w(x, t)\ 

and 



R'' 



\{E(w{x,s)) ■ (u{x,s)^w(x,s)))\e 2 Ss\\ s M\\Ur)dr dx 



ds 



lim 

v->0 



R ' 



\(E(w{x, s)) ■ (u v {x, s) - w(x, s)))\dx e 2 Si WS^U^dr 



ds . 



Consequently, one has: 



\u(x,t) - w(x,t)\ 2 < ( / |uo(x) ~w(x,0)\ 2 dx )e 2 So ISMIUW* 



+ 2 



K^Mx, s)) • (u(x, s) - w{x, s)))dx\e 2 SI ISMIU(r)rfr 



R'' 



(60) 



ds . 



Remark 4.3 1. The most important remark following this proof is the fact 
that in the presence of a smooth solution of Euler equations any weak limit 
will coincide with this smooth solution and since the smooth solution conserves 
energy one has: 



lim \v \ \ \\Iu v (x, s)\dxds 
v Jo Jr" 



<l f \u (x)\ 2 dx- limi f \u„(x,t)\ 2 dx (61) 
<\ I \Mx)\ 2 dx-\l \u(x,t)\ 2 dx = 



and therefore there is no anomalous dissipation of energy ( the energy dissipation 
rate which appears in Kolmogorov theory of turbulence is 0.) 
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2. In fact, as said above, regularity is known to propagate for smooth solution of 
the 2d Euler equations, therefore, in 2d with smooth initial data the energy will 
always be conserved and there will be no anomalous dissipation of energy (or 
enstrophy). The appearance of singularities, after some (long?) time, for the 
3d Euler equations with smooth initial data is still a challenging open problem. 
Therefore, it may happen (even with smooth initial data) that, after some long 
time, anomalous dissipation of energy appears in the limit, as v — > 0. Thus, 
anomalous dissipation would be a criterion for blow up of the solutions of 3d 
Euler equations. 

3. Moreover, returning to the contributions of De Lellis and Szekelyhidi, one 
should keep in mind that there exists in 2d or 3d an infinite set of initial data, 
uq € L 2 , such that for each Uq there is an infinite set of corresponding dissipative 
solutions for Euler equations. In this situation the above consideration does not 
apply. For such initial data there may be anomalous dissipation of energy. In 
this case, since the fact of being dissipative does not imply the uniqueness of 
the solution, one should consider other types of uniqueness criteria, or selection 
mechanism for uniqueness. A natural candidate for unique solution of Euler 
equations would be a viscosity limit solution, namely the limit solution, u, of a 
family of solutions, u v , of Navier-Stokes equations, when the viscosity v — > 0. 
This assertion is not proven with full generality, but some examples built with 
the shear flow, cf JIOf . indicate that this should be the case. 

Remark 4.4 Extension of the Proposition ^. 3\ assuming more regularity on the 
initial data, can be found in J67jj (see also JEEl). Assuming that Euler equations 
have a smooth solution on the time interval [0,T], it is shown that the same is 
true for the solution of Navier-Stokes equations with the same initial data on 
the same time interval for v small enough. Moreover, the rates of convergence 
in higher norms are estimated (see also \59f for similar results concerning the 
Euler— a regularization model.) 

4.4 Viscosity limit with boundary effects 

Unlike the previous section the emphasis here is put on the viscosity limit, in 
both the 2d and the 3d cases, in the presence of physical boundaries. In fact, in 
as far as our understanding of this problem goes, there is not much difference 
between the 2d and the 3d cases. Thus, it will be assumed that the solution 
of Euler equations (with given fixed and smooth initial data uq) is smooth for 
t € [0, T]. In the 2d case the global regularity of Euler equations, i.e. for all 
T > 0, was proven by Wolibncr in 1933 [88] (and independently at the same time 
by E. Holder [48]), and it will be assumed in the present discussion for the 3d 
case as well. As said above, the introduction of Navier-Stokes equations allows 
to take into account additional boundary effects (in particular the production 
and shedding of vorticity near the boundary). One may link these effects with 
anomalous dissipation of energy, as v — > 0, and conclude that turbulent regimes 
which are observed in the fluid flow, even far from the boundary, where one 
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could talk of homogenous and isotropic turbulence, is not a smooth solution of 
Euler equations. 

In fact, the number of boundary conditions required for Navicr-Stokes equa- 
tions is 3, and this number is reduced to 1 in the case of Euler equations, i.e., 
u ■ ft = 0, generating a boundary layer near the boundary . However, since the 
problem is nonlinear, this boundary layer may not be confined to the boundary 
and instead it may propagate, by means of the nonlinear advection term, inside 
the bulk of the flow. 

Here also dissipative solutions to the Euler equations will be used; but an 
important modification will be introduced in their definition. In the present case, 
smooth solutions of Euler equations may have nonzero tangential component 
at the physical boundary, thus comparing them to test functions that have 
compact support in the domain (i.e. without interaction with the boundary) is 
not enough to characterize the limit of Navier- Stokes solutions, as the viscosity 
v — > 0. The modified definition of dissipative solutions for the Euler equations 
is: 

Definition 4.1 Consider the Euler equations in a domain CI, with nonempty 
boundary dQ 7^ 0. A vector field u is called a dissipative solution up to the 
boundary, for the Euler equations, if for every smooth divergence- free test vector 
field, w, which is tangent to the boundary dfl (i.e., w ■ ft = on dCl), and which 
has bounded support in Q, one has the estimate: 



where E{w) is the residual as in equation i26}) . 

This is the same criterion as in (|27]l. but with a larger class of test functions, w. 
In particular, if there exists a smooth solution u(x,t) with initial data uq(x), 
any dissipative solution up to the boundary with the same initial data does 
coincide with this smooth solution. Such a result would not be valid if only test 
functions with compact support in f2 were used in the criterion (|27[) . 

Theorem 4.1 In dimensions 2 or 3, for any v dependent family of Leray-Hopf 
weak solutions, on the time interval [0, T], of the Navier-Stokes equations 




(62) 



d t U V + Uy ■ Vli; 



l/AU V + Vf>; 



in a domain Q, with boundary conditions 



u y ■ ft = , v(pfiUy + C(x)u v ) T + \(v, x)u v = on dfl , 
and X(u, x)>0, C{x) e C(R n ; R n ) ; 



(63) 
(64) 
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and with fixed initial data uq, one has the following assertions: 

1. Ifu v converges weak—*, in L°°((0, T); L 2 {0)) ; to a function u(x,t) , which 
is, for all t G (0, T), a smooth solution of the Euler equations, then u v converges 
strongly to u, and it does not dissipate energy in the limit as v — > 0; namely, 

eO,T)=lim / {v [ \Vu v \ 2 dx + [ \(v)\u v (x, t)\ 2 da)dt = . (65) 
Jn Jdn 

2. If (y@Fpjr) T — > 0, as v — y 0, even in a very weak sense (for instance, in the 
sense of distribution in T>' (dfl x (0, T)) ), then any weak limit, u, of the sequence 
u v is a dissipative solution up to the boundary of the Euler equations. 

3. (Kato criterion) If the sequence u v satisfies the estimate: 

lim [ v [ \Vu v (x,t)\ 2 dxdt = 0, (66) 

then the hypothesis of the point 2 above is satisfied, and therefore any weak limit 
of the sequence u u is a dissipative solution up to the boundary of Euler equations. 

Proof Below, we will sketch some of the main ideas of the proof. 

Equation (|6"5)) represents the dissipation of energy over the time [0,T], the 
term J^{y J n \ Vu L ,\ 2 dx)dt represents the dissipation of energy due to the fluid 

viscosity, and (f 9n \(i>)\u u (x,t)\ 2 do-)dt represents the dissipation of energy 
due to the friction of the fluid flow with the boundary. In the special case of the 
Dirichlet boundary condition for the Navier- Stokes equations this second term 
does not appear because the velocity field, u v , is zero on the boundary. 
Point 1 follows from the energy estimate 

^ f \u v {x,t)\ 2 dx+v I [ \Vu u \ 2 dx+ [ \(x)\u u (x,t)\ 2 da 



a Jo Jn Jon 

< [ \u \ 2 dx + o(l), 
Jn 

where o(l) stands for a term that tends to zero, as v — > 0. The above energy 
inequality is obtainedin a similar way as for the Leray-Hopf weak solutions of 
the Navier-Stokes system; however, this time using the boundary conditions 
(|63 ]) - (|64| . For this reason we also call weak solutions of Navier-Stokes equations 
that satisfy the above energy inequality Leray-Hopf weak solutions. 

Observe, thanks to the weak—* convergence in L°°((0, T); L 2 (yi)) of u v to a 
smooth solution u of Euler equations, that one has: 

- I \u {x))\ 2 dx = - I \u(x,t)\ 2 dx = ]- t \u{x,t)\ 2 dx 



i r (67) 

< lim — / \u v (x,t)\ dx . 
u^o 2 J n 



Concerning point 2, the main difference between the problem without phys- 
ical boundary and the current problem with boundary comes from the, only 
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uncontrolled, extra term in the computation of 

rT 



o Jn 



(Ait„ • {u v — w))dxdt . 



which, thanks to the boundary conditions w ■ n — and (163p ~ (|64[) . gives a 
boundary term of the form 

rT 



((^p-) T ' w T )dadt. 
o JdQ. an 



(68) 



Thus to establish point 2 one needs to show that the above term goes to zero, 
as v — >• 0. By virtue of the density of the set of smooth test functions, to, in 
Definition 14 . 1 1 (sec estimate ([62]) ). it is enough to consider in the criterion stated 
in point 2, and to test it for, C°° test functions. 

Next, we show point 3. For any smooth w{x,t) tangent to dil for all t € 
(0, T], one constructs (following Kato [51]) a family w v (s,T,t) (here r and s are 
the local tangential and normal coordinates at dfl, respectively), for small v > 0, 
of smooth divergence-free test functions, which arc tangent to the boundary, <9f2, 
and which also satisfy the following properties: For every v > 0, w v coincides 
with w on dtt x (0,T), and its support is contained in 

Sl v x (0, T) = [i£!l: d(x, dn) < v} x (0, T) . 

Moreover, w v satisfied the estimates 

C 



\d t w v \ L °° < C , iV-rttvlioo < C, \d s w u \ L <^ < 



(69) 



Next, we show that (|66f implies \\m v ^{{v^s-) T , uv)z, 2 ((9nx(o .t)) = 0. Ob- 
server that ((^^)r,W T ) L 2( anx ( 0!T)) = {{v^)r,{w v )r)L3{dSlx(fi,T))i and henCe 

this expression can be deduced from Navier-Stokes equations, subject to the 
boundary conditions (|63| -([64 |) . according to the formula: 







{0,w v )dt 



{u u ,d t w u )dt 



{{d t u v + V ■ (u v ® u v ) — v>Au v + Vp„), to y )df 



((mk • V)w„,u v )dt 



(70) 



Jo 

which implies the estimate: 



(vdnu u , uv)l 2 (3Qx(o,t)) 



(K • V)uv, u v )dt 



o(l) 



(71) 



Eventually, using (j6"9")l . the control of the energy and the Poincare inequality, 
one obtains 



((u v ■ V)to„, u v )dt 



< C 



J Q 



(72) 
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Therefore, thanks to ([7T|) and (|T2")) one concludes 

^o^¥' r,WT ^ J(flnx(ll ' T)) = °' 
This completes the proof of point 3. 

Remark 4.5 The boundary conditions l63\ )- [64\ ) together with the hypothesis 
lim„_5.o A(z/) = imply that any weak limit of the family u v , of Leray-Hopf 
weak solutions of the Navier-Stokes equations, is a dissipative solution up to the 
boundary of Euler equations. The hypothesis of point 2 involves only v times 
the tangential component of the normal derivative of the velocity. This point 
of view can be combined with the hypothesis of point 3, and this was done by 
Wang '851. The condition proposed by Wang (which at the end of the day will 
be equivalent to all the others) involves only the tangential component of the 
normal derivative of the velocity, but in a layer slightly thicker than the Kato 
layer, which is the well-known thickness of viscous sublayer as predicted by von 
Karman (1930) logarithmic profile. In situation when the convergence to the 
regular solution of Euler equations is insured it is possible to give many details 
of the rate of this convergence, the nature of the boundary layer and so on. 
There is a huge literature along these lines, to cite few see, for example, \13f . 
\mi , fgffi and references therein. 

Remark 4.6 When A(V) does not go to zero, in particular, in the most extreme 
case (the Dirichlet boundary conditions) when X(v) = oo, there is no information 
on the nature of the limit, at this level of generality. 

However, one can underline the interpretation of the above theorem by stating 
that any of the following facts: 

• There is no production of vorticity at the boundary. 

• There is no anomalous dissipation of energy in a layer of size v , at the 
boundary, and, of course, no anomalous dissipation of energy in the whole fluid 
flow domain. 

• There is conservation of energy as v — > , 

implies that the zero viscosity limit does not exhibit any turbulent behavior. 
This is meant in the sense that the limit solution of Navier-Stokes equations, as 
v — > 0, is as smooth as the corresponding solution of Euler equations, and that 
no nontrivial Reynolds stresses tensor (in the sense of coarse graining) appears 
at the limit. 

On the other hand, if one of the above statement is not true, due to boundary 
effects, one cannot say that the limit is a dissipative solution up to the boundary 
of Euler equations in contrast with the case without boundaries. 

Nontrivial Reynolds stresses tensor (in the sense of coarse graining) may 
occur either when the weak limit, u, as v — > 0, is not a weak solution of Euler 
equations, or whenu is a weak solution of Euler equations that is obtained as the 
limit of Navier-Stokes equations with energy decay, i.e., u is a weak admissible 
solution of Euler equations. However, even so, some concentration of oscilla- 
tions (turbulence) near the boundary may prevent it to be a dissipative solution 
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up to the boundary of Euler equations. This pathologic type of weak solutions of 
Euler equations can in fact be constructed following J 11}/ . 



4.5 More about boundary effects and Prandtl equations 

In this section we follow very closely our presentation (cf. [8]) concerning the 
Prandtl equations. What essentially prevents the existence of more definite 
results, in the general case, is that what would be a boundary layer, in the 
zero viscosity limit, may not remain near the boundary and will actually detach 
from the boundary due to its rolling up into vortices which will, by the advection 
term, penetrate the bulk of the fluid flow forming the wake of the solid body. 

However, this fact does not prevent more specific theoretical or practical 
studies. If one assumes, as a hypothesis, that the variation of the tangential 
component, and the size of the normal component, of the velocity of the fluid, 
near the boundary, are of the order of e~ 1 and e, respectively, (with e = 
then by a change of scale near this boundary one may use the Prandtl equations 
(cf. |73| and for more recent discussion concerning the relation between the issue 
of viscosity limit and the Prandlt equations [86]) We recall below the planar 
Prandtl equations for the case of a flat boundary: start from the Navier-Stokcs 
equations in the half-plane, {(xi,x 2 ) <E R 2 : x 2 > 0}, with the no slip boundary 
condition, u v (xi,0,t) = 0, 

d t ul - v^u\ + u\d x ,u\ + u v 2 d X2 u\ + d Xl p v = , (73) 

d t u v 2 - i/Ai£ + u\d Xl u v 2 + u v 2 d X2 u v 2 + d X2 p v = , (74) 

d Xl u\ + d X2 u v 2 =Q, (75) 

<(xi,0) = tiafci.O) = for i^R, (76) 

then introduce the scale e = yjv and the change of variable which corresponds 
to a boundary layer in a parabolic problem, taking into account that the normal 
component of the velocity on the boundary, U2, remains 0: 

X 1 = x u X 2 =^, (77) 
Ul(xi,X 2 ) = ui(xi,x 2 ),u 2 (x 1 ,X 2 ) = eu 2 (x 1 ,x 2 ) . (78) 

Going back to the original notation {x\, x 2 ) and letting e go to zero, one obtains 
formally the system: 

u 1 (a;i,0 ) *)-f7i(xi,i)=0, (79) 

d X2 p(xi,x 2 ) = p(x 1 ,x 2 ,i) = P{xi,t) (80) 

d t ui - dl 2 ui + uid Xl U! + u 2 d X2 u x = -d Xl P(xi,t) , (81) 

d Xl ux + 3 X2 u 2 = , (82) 

ui(a;i,0,t) = u 2 (xi,0,t) = for neR, (83) 

lim ui(xi,x 2 ) = 0. (84) 

X2— >00 
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Equations ([75)1 and describe the fact that the boundary layer solution 
will match a solution of incompressible Euler equations. U\{xi,t) and P{x\,t) 
are the tangential velocity and the pressure, respectively, on the boundary, 
X2 = 0, of the solution to the Euler equations with the same initial data. Finally, 
the solution is written as a sum of regular functions: the boundary layer term, 
as computed above, and the solution of Euler equations, denoted u- ln t (xi, x 2 , t) 



V(*1,*2,*A = ( uKx^^t) N 

u^{x 1 ,x 2l t)J \yj)u v 2 {x u ^,t) 



+ u iat (xi,X2,t) . (85) 



Remark 4.7 The validity of the Prandtl approximation is consistent with the 
Koto's criterion (point 3 of Theorem ^. J[ ) because with &85]) one has: 

|V A u„(x, s)\ 2 dxds < Cy/v . (86) 

Jnr\{d{x,dQ,)<cv} 

Remark 4.8 At least one example showing that the Prandtl expansion cannot 
be always valid has been constructed by Grenier f^7| /. Solutions are considered 
in the domain 

R Cl /ZxR+ . 

Grenier starts with a solution u v re ^ of the pressure-less Navier-Stokes equations 
given by: 

< e/ = K e/ (t,4=).0), (87) 

d t U re f - OyUref = , (88) 

where Y = . Using a convenient and explicit choice of the function u re f , 
with some sharp results on instabilities, a solution of Euler equations of the form 

u = u ref + Sv + 0(5 2 e 2Xt ) for < t < — — - (89) 

A|logi5| 

is constructed. It is then shown that the vorticity generated at the boundary, by 
the solution of Navier-Stokes equations (with the same initial data), is too strong 
to allow the convergence of the Prandtl expansion. One should observe that, once 
again, this is an example with infinite energy. It would be interesting to see if 
such example could be modified to enter the class of finite energy solutions and 
analyze how it would violate the Koto's criterion (66\) . 

It is important to observe that the mathematical properties of Prandtl equations 
exhibit the pathology of a situation based on the hypothesis that the boundary 
layer does not detach. First, one can prove the following 

Proposition 4.1 Let T > a finite positive time, and let (U(x,t), P(x,t)) € 
C 2+Q (Rt x (R X1 x Ri)) be a smooth solution of the 2d Euler equation satisfying 
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for t = the compatibility condition £/i(xi,0, t) = C/ 2 (xi,0,i) = (only the 
boundary condition [/ 2 (xi,0,i) = is preserved by the Euler's dynamics). Then 
the following facts are equivalent: 

1. with initial data u = 0, the boundary condition f2i(xi, 0, f) = Ui(x\,0,t) 
in J7ff| ) and the right-hand side in i8U\) given by P(xi,t) = P(xi,0,t) , Prandtl 
equations have a smooth solution for < t < T . 

2. The solution u v {x, t), with initial data u„(x, 0) = U (x, 0) and with no-slip 
boundary condition, at the boundary x\ = 0, converges in C 2+a to the solution 
of Euler equations, as v — > 0. 

The fact that the hypotheses 1. (and 2.), above, may be violated for some t is 
related to the appearance of a detachment zone and the generation of vorticity 
at the boundary (which may eventually generate turbulence in the bulk of the 
fluid flow) . This is well illustrated in the analysis of Prandtl equations written 
in the following simplified form: 

d t ui - dl 2 ui + uid Xl ui + u 2 d X2 ui = , d Xl u\ + d X2 u 2 = , (90) 
ui(xi,0) = £t 2 (xi,0) = for xi G R, (91) 
lim ui(xi,X2)=0, (92) 

X2— >O0 

ui(xi,x 2 ,0) = uo(xi,x 2 ) . (93) 

Regularity in the absence of detachment corresponds to a theorem of Oleinik 
[70] . She proved that global smooth solution does exist if the initial profile is 
monotonic, i.e. 

m(x,0) - (ui(xi,x 2 ),0) ,d X2 ux(x 1 ,x 2 ) ^ 0. (94) 

On the other hand, initial conditions with "recirculation properties" leading to a 
finite time blow up have been constructed by E and Enquist [35] . An interesting 
aspect of these examples is that the blow up generally does not occur exactly 
on the boundary, but inside the domain supporting the hypothesis that this 
blow up is due to the roll up of the boundary layer which increases the velocity 
gradients. 

The above pathology appears in the fact that the Prandtl equations are 
highly unstable. This comes from the determination 112 in term of u\ by the 
equation: 

d Xl ui + d X2 u 2 = . (95) 

With this observation more recent and systematic studies of these instabilities 
can be found in [43] and |44j . 

Therefore it is only with analytic initial data (in fact analyticity with respect 
to the tangential variable is enough) that one can obtain (using an abstract 
version of the Cauchy-Kowalewsky theorem) the existence of a smooth solution 
of Prandtl equations for a finite time and the convergence to the solution of Euler 
equation during the same time (cf. Asano, Caflisch-Sammartino and Cannone- 
Lombardo-Sammartino [2], [75] . [61] . and more recently with less stringent 
hypothesis on the x 2 asymptotic behavior in |53j.) 
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5 Conclusion 



This contribution is organized around the issues of the Cauchy problem, and 
the relevant boundary effects, in equations of motion of incompressible fluids. 
It presents some mathematical tools and results (theorems) that have been de- 
veloped over the last, and the present, century for investigating certain basic 
properties of solutions of equations of fluid flows. There are many different 
mathematical approaches, that complement each other, to address these chal- 
lenging problems; however, we chose to focus, in this contribution, only on 
certain mathematical aspects that are concerned with: 

• The relation between the different equations of fluid motion in terms of 
well identified scaling parameters. 

• The existing results (far from being complete at the time of this writing) 
about classical solutions of these equations. 

• The notion of weak solutions (as used by Leray in his thesis or in the 
sense of distributions, as introduced by Gelfand and Schwartz). In particular, 
we recall that the notion of weak solutions (even with conservation or decay of 
energy) is not precise enough to characterize the convenient (physically relevant) 
solutions of the Cauchy problem for Euler equations. This is in spite the fact 
that the relation between the regularity of solutions of Euler equations and the 
non conservation of energy is related to the Onsager conjecture in turbulence 
theory. 

• The analysis of the zero viscosity limit for solutions of Navier-Stokes equa- 
tions, with or without boundary effects. In the absence of boundary effects 
this convergence is a direct consequence of the existence (or nonexistence) of 
smooth solutions, with the same initial data, for the Euler equations; while in 
the presence of boundary effects the situation is completely different and the 
non convergence (following a theorem of T. Kato) is related to the anomalous 
dissipation of energy, in the limit as v — > 0. This anomalous dissipation of en- 
ergy corresponds, in the Kolmogorov theory of turbulence, to the non- vanishing 
limit of the positive mean rate of dissipation of energy, e, as v — > 0. 

• Lastly, the Prandtl equations are shortly discussed. We recall that these 
equations are in general ill-posed and that in the case when they are well-posed, 
which is the "laminar regime" , the Kato criterion (of the above mentioned Kato 
theorem) is obviously satisfied, and as a result there is no anomalous dissipation 
of energy. Beside their importance (in situations where they are valid) for 
engineering sciences applications, the Prandtl equations became basic models for 
a systematic mathematical analysis of boundary layer in many other situations. 
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